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TABLES 


I.  INTRODUCTION 


In  a previous  report,  band  model  methods  were  used  to  compute  the 
effective  transmittance  of  the  atmosphere  for  radiation  emitted  by  a high- 
temperature  gaseous  H2O/CO2  source  with  characteristics  representative  of 
a missile  plume  at  an  attitude  of  20  km.  The  band  model  formulation  used 
allowed  for  the  high  degree  of  correlation  that  exists  between  the  line  positions 
of  the  source  emission  spectrum  and  the  atmospheric  H20/C02  absorption 
spectrum  by  treating  the  entire  line  of  sight  through  the  atmosphere  and  the 
hot-gas  emission  source  as  a single,  highly  inhomogeneous  optical  path. 
Inhomogeneities  along  the  path  were  accounted  for  with  either  the  Curtis  - 
Godson  (CG)  or  the  Lindquist-Simmons  (LS)  approximation. 

The  band  model  formulation  was,  however,  limited  by  the  assumption 
that  the  individual  spectral  lines  within  the  band  could  be  described  by  a 
pressure -broadened  Lorentz  profile.  The  effect  of  this  limitation  became 
apparent  in  subsequent  calculations  made  for  hot  gas  sources  that  represented 
missile  plumes  in  the  sea  level  to  60  km  altitude  range.  As  the  target  altitude 
increased  above  about  20  km,  the  computed  effective  transmittance  to  space 
approached  unity  much  more  slowly  than  the  transmittance  predicted  by  con- 
ventional atmospheric  band  models.  Although  a substantial  part  of  the 
apparent  enhancement  of  the  atmospheric  absorption  is  caused  by  line  corre- 
lation effects,  much  of  the  effect  was  also  attributed  to  the  inappropriateness  of 
using  a pressure-broadened  Lorentz  line  shape  for  such  high  altitudes. 

Typical  line  widths  for  pressure  broadening  in  atmospheric  applications 
are  of  the  order  y^  ~ (0.  1 cm  /atm)  p,  where  p is  the  total  gas  pressure.  At 
an  altitude  of  25  km,  the  pressure  is  p ~0. 03  atm,  and  y ^ ~ 0.003  cm  ^ . For 
higher  altitudes,  both  p and  y^  are  smaller.  The  mechanism  of  Doppler 


S.  J.  Young,  Band  Model  Calculations  of  Atmospheric  Transmittance  for  Hot 
Gas  Line  Emission  Sources,  TR-007 5(5647 ) - 1 , The  Aerospace  Corporation, 
El  Segundo,  Calif.  (31  December  1974). 
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broadening,  on  the  other  hand,  is  pressure  independent  and  depends  only  on 
the  local  temperature,  the  mass  oi  the  absorbing  or  emitting  molecule,  and  the 
wavelength  of  the  spectral  line.  For  CO?  in  the  2. 7 -pm  region,  for  example, 
Yd~  0.003  cm  throughout  the  whole  atmosphere.  Thus,  above~25km,  the 
use  of  a pres  sure  - broadened  Lorentz  line  will  give  a smaller  line  width  (and 
different  shape)  than  actually  exists.  Since  any  atmospheric  slant  path  that 
extends  to  space  has  some  portion  above  25  km,  Doppler  broadening  should  be 
accounted  for.  Because  of  the  exponential  decrease  of  pressure  and  absorber 
concentration  with  altitude,  most  of  the  radiative  transfer  effects  for  an 
atmospheric  path  will  occur  along  the  low-altitude  end.  If  enough  of  the  path 
is  below  25  km,  then  the  use  of  a pressure -broadened  Dorentz  line  for  the 
whole  path  may  be  justified.  It  is  only  for  paths  that  lie  totally  above  ~2  5 km 
that  Doppler  broadening  must  be  accounted  for. 

Both  pressure  broadening  and  Doppler  broadening  can  be  accounted  for 
by  using  the  Voigt  line  shape.  This  line  shape  is  a convolution  of  the  dispersion 
(Lorentz)  and  Gaussian  (Doppler)  line  shapes  and  reduces  to  either  of  these 
component  line  shapes  in  the  appropriate  temperature  and  pressure  limits. 

This  report  is  concerned  with  the  use  of  the  Voigt  line  shape  in  band  model 
calculations  of  the  effective  transmittance  of  the  atmosphere  to  hot  gas  line 
emission  sources. 

A review  of  the  band  model  formulation  used  in  the  computation  of  the 
effective  transmittance  and  the  approximation  used  to  incorporate  the  Voigt 
line  shape  are  presented  in  Section  II.  In  Section  III,  results  obtained  with  the 
CG  and  LS  approximations  are  presented,  and  the  limitations  of  the  approxi- 
mations are  discussed.  In  Section  IV,  an  improved  inhomogeneity  approxi- 
mation is  formulated,  and  the  results  of  the  formulation  are  presented. 


Collisiona)  narrowing  effects  on  the  Doppler  line  shape  and  width  are  not 
considerei  here.  The  whole  concept  of  non-Gaussian  shapes  and  pressure 
dependence  of  widths  of  Doppler  lines  has  only  recently^  been  introduced  into 
the  field  of  atmospheric  transmittance  and  absorptance  computations. 

R.  L.  Armstrong,  ’'Effect  of  Collision  Narrowing  on  Atmospheric 
Transmittance,"  Appl,  Opt.  14,  56  (197  5). 


-6- 


-r^'- — 


temperature  T(s),  and  r (v  , s)  is  the  mean  transmittance  for  av  over  the  path 
between  s = 0 and  the  general  position  s.  If  more  than  one  absorbing  molecular 
species  exists  in  the  path,  the  total  transmittance  is  assumed  to  be  the  product 
of  the  transmittances  computed  for  the  individual  species;  i.e.,  it  is  assumed 
that  no  line  correlation  exists  between  the  absorption  spectra  of  different 
species.  For  the  statistical  band  model  formulation  in  which  the  number  of 
lines  N in  Av  approaches  infinity,  and  for  which  the  assumption  is  made  that 
the  absorption  caused  by  each  line  whose  center  occurs  in  Av  is  confined  to 
Av  , the  mean  transmittance  for  a single  species  is 


r(s)  = exp 




3 

S.  J.  Young,  Band  Model  Formulation  for  Inhomogeneous  Optical  Paths, 
TR-0075(5647)-4,  The  Aerospace  Corporation,  El  Segundo,  Calif. 

(19  December  1974);  also,  J.  Quant.  Spectr.  Radiative  Transfer  15,  483  (197^). 

4 

S.  J.  Young,  Addendum  to;  Band  Model  Formulation  for  Inhomogeneous 
Optical  Paths,  TR-O076(697O)- 3,  The  Aerospace  Corporation,  El  Segundo, 

Calif.  (31  July  1975). 
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where  W is  the  mean  equivalent  width  ol'  all  the  lines  in  Av  , and  6 is  the  mean 
line  spacing  parameter  6 Av  /N.  The  functional  dependence  on  v and  species 
is  suppressed  in  Eq.  (2)  and  from  here  on.  Evaluation  of  the  integral  of  Eq . (1) 
requires  the  transmittance  derivative,  which,  from  Eq.  (2),  is 


dT  (s) 
ds 


t(s)  - --W(s) 
6 ds 


(3) 


The  mean  equivalent  width  derivative  can  be  written  in  the  form 

i = c(s)p(s)k  (s)y[x(s),  p(s)l 


(4) 


where  c(s)  and  p(s)  are  the  mole  fraction  of  absorber  and  total  pressur<  distri- 
bution, respectively,  along  the  path,  and  k(s)  is  the  mean  absorption  coefficient 
band  model  parameter  for  the  interval  Av  and  temperature  T(s).  The  argu- 
ments of  the  function  y(x,  p)  are 


x(s) 


k (s)u{s) 

Ms) 


and 


p(s)  = 


P (s) 

{3e(s) 


(6) 


whe  re 


ke,5) 


1 

u(s) 


c(s/)p(s/)k(s')ds/ 


(7) 
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and 


{3  (s)  = ^ f c(s/)p(s/)k(s/)p  (s')ds' 

e u(s)Ee(s)  JQ 


(8) 


are  path-averaged  values  of  k and  (3,  respectively.  The  parameter  x(s)  is  a 
measure  of  optical  depth  between  s = 0 and  s,  and  p(s)  is  a measure  of  the  local 
inhomogeneity  at  s.  The  mean  line  width  to  spacing  parameter  (3  is  defined  by 


\L(s) 


(3l(s)  = 2tt  g (g) 


(°>a) 


for  Lorentz  lines  and  by 


pd(s)  = 


1 i,  Vs> 

In  2 6 (s) 

e 


(9b) 


for  Doppler  lines.  In  Eq.  (9a),  "y  s ) is  the  mean  Lorentz  line  width  for  all 
the  lines  in  Av  and  for  the  temperature,  pressure,  and  species  concentration 
conditions  at  s.  In  Eq.  (9b),  -y  ^( s ) is  the  mean  Doppler  line  width  for  Av  and 
the  temperature  T(s).  In  either  equation,  6g(s)  is  the  effective  mean  line  spacing 
band  model  parameter  for  Av  and  temperature  T(s).  The  optical  depth  is 


u(s)  = / c(s/)p(s/)ds/ 


(10) 


The  functional  form  of  y(x,  p)  depends  on  the  assumed  line  shape,  the 
assumed  distribution  of  line  strengths  in  Av , and  the  approximation  used  to 
account  for  inhomogeneities  along  the  path.  Only  the  exponential -tailed  inverse 
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line  strength  distribution  is  considered  here.  With  this  limitation,  the  functional 
form  for  the  Lorentz  line  shape  within  the  CG  approximation  is 


y(x,  p)  = (2 


+ (1  - p) 


f(x) 


X 


where 


f(x)  = ^ (^1  + 2ttx  - 1 ) 


(in 


(12) 


is  the  equivalent  width  curve  of  growth  function.  In  the  LS  approximation,  the 
result  for  the  Lorentz  line  shape  is 

V(K.D)  ^ 2P(‘  tpZ)sg4JlX.  (13) 

yi  + 2ttx  (p  + yi  + 2ttx) 

The  solutions  of  y(x,  p)  for  Doppler  lines  are  not  expressible  in  closed  form,  but 
are  conveniently  obtained  by  tabular  interpolation  and  approximate  solutions. 

For  the  LS  approximation,  the  data  of  Table  1 can  be  used  to  obtain  y(x,  p)  bv 
interpolation  for  most  practical  applications.  For  small  x,  the  series  solution 


°o  n 

y(x,  p)  = £ 

n = 0 vl  + np^ 


The  additional  factors  of  tt  that  appear  in  Eq.  (13)  [and  later  in  Eq.  (17)],  but 
not  in  Eq.  (9)  or  Eq.  (6)  of  Ref.  4,  are  included  in  order  to  compensate  for  the 
form  of  the  band  model  parameter  6e  used  with  the  band  model  (Section  III). 
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is  convenient,  whereas,  for  p < 0.5,  the  approximation 


y(x,  p) 


2 

j(l  + TTX)[l  - TTX(P2  - l)]j1/2 


(1?) 


is  accurate  to  better  than  1 percent.  In  the  CG  approximation,  yfx,  p)  is  defined 
by 


y(x,p)  = (2  - p )^£>  + (1  - p)^T  (16) 

whe  re 

H(x)  =?r2/”ln 


1 + TTXO 


-z 


dz 


(17) 


is  the  equivalent  width  curve  of  growth  function  for  Doppler  lines  and  the 
exponential- tailed  inverse  line  strength  distribution.  For  small  x 


HW  = i 2 


(-TTX) 


n+1 


n = 0 (n  + 1 ) 


3/2 


(18) 


whereas,  for  large  x 

1 , *2/8  , 7 tt4 / 640 

2 T 4 v ' 

(In  ttx)  (In  ttx) 

For  intermediate*  x,  H(x)  can  be  interpolated  from  the  data  of  Table  2.  The 
derivative  dH(x)/dx  may  be  determined  by  differentiating  Eqs.  (18)  and  (19)  or 
by  performing  a numerical  differentiation  of  the  data  in  Table  2.  In  practice, 


H(x)  ^ — —j = (In  ttx)3^2 

3*  16 
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Table  2.  Curve -of -Growth  Function  for  Doppler  Line 
Shape  and  Exponential-Tailed  Inverse  Line 
Strength  Distribution 


TTX 

ttH(x) 

0.  1 

0.09665 

0.2 

0.  18722 

0.4 

0.35325 

0.6 

0. 50289 

0.8 

0.63940 

1.0 

0.76515 

2.0 

1.28138 

4.0 

2.00567 

6.0 

2. 52794 

8.0 

2.94136 

10 

3.28568 

20 

4.47964 

40 

5.84661 

60 

6.71916 

80 

7.36866 

100 

7.88910 

200 

9.59377 

400 

11.4233 

600 

12. 5478 

800 

13.3689 

1000 

14.0187 
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it  is  more  convenient  to  obtain  dH(x)/dx  for  intermediate  x from  the  data  of 

Table  1 by  interpolating  on  the  p = 1 column.  This  is  possible  because,  for 

P = !’  i,e*'  for  a homogeneous  path,  the  CG  and  LS  solutions  for  y(x,  1)  are 
identical,  and  for  either 


Y(x,  1) 


dH(x) 

dx 


(20) 


The  equivalent  width  factor  of  Eq . (3)  is  given  by  Eqs.  (4)  through  (20). 
The  transmittance  factor  of  Eq.  (3)  is  given  by  Eq.  (2),  where,  for  the  CG 
approximation 


W(s) 

5 


Pe(s)f  (x(s)] 


(21a) 


for  Lorentz  lines  and 


W(s) 

5 


Pc(s)Hfx(s)] 


(21b) 


for  Doppler  lines.  In  the  LS  approximation, 
derivative  function 


W(s)/6  must  be  integrated  from  the 


W(s)  _ 
5 " 


1 dW(s')  . , 
6 ds  ds 


(22) 


for  either  line  shape. 

The  formal  representation  of  the  statistical  band  model  method  for  com- 
puting the  radiance  of  inhomogeneous  optical  paths  is  given  by  Eqs.  (1)  through 
(22)  for  Lorentz  or  Doppler  line  shape  in  either  the  CG  or  LS  approximation. 
Some  comments  are  now  given  concerning  the  application  of  the  model  to  the 
present  problem,  where  the  optical  path  is  composed  of  two  distinct  component 
paths.  Let  the  path  position  s - 0 represent  the  point  where  the  spectral 
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radiance  is  being  measured.  (For  a space  sensor  geometry,  s = 0 is  the 
altitude  above  which  the  atmosphere  is  assumed  to  be  nonabsorbing.  ) Let 
s = Sq  be  the  boundary  between  the  atmosphere  and  the  near  (to  the  sensor) 
boundary  of  the  hot  gas  source,  and  s = S represent  the  far  boundary  of  the 
source.  If  Eq.  (1)  is  written  as  the  sum  of  integrals  from  s = 0 to  sQ  and  from 
Sg  to  S,  then  the  first  term  represents  the  atmospheric  contribution  to  the 
measured  radiance,  and  the  second  term  represents  the  hot  gas  source  con- 
tribution to  the  radiance  with  allowance  made  for  absorption  by  the  intervening 
atmosphere.  All  results  presented  here  are  for  the  2.7-|j.m  region,  and,  for 
the  target  plume  temperatures  used,  the  atmospheric  contribution  to  the  total 
radiance  can  be  neglected.  The  second  term 


(23) 


also  represents  the  unattenuated  radiance  of  the  source  if  the  intervening 
atmosphere  is  nonabsorbing.  The  effective  transmittance  of  the  atmosphere 
is  computed  by  taking  the  ratio  of  the  solution  of  Eq.  (23)  for  the  absorbing 
atmosphere  case  to  the  solution  for  the  nonabsorbing  case,  and  is  designated 
Te*  Because  this  method  of  computing  transmittance  implicity  takes  line 
correlation  into  account,  will  always  be  smaller  than  the  transmittance  t 
computed  from  Eqs.  (2)  and  (21),  or  (22),  with  s = sQ. 

The  two  conditions  for  which  Eq.  (23)  must  be  evaluated  are  conveniently 

accommodated  by  writing  u(s),  k (s),  and  (3  (s)  for  s > sn  as 

0 0 * 0 


u(s)  = u( Sq ) + 


c(s  )p(s/)ds/ 


(24) 


ke<B> 


1 

u(s) 


U(VW + 


c(s/)p(s/)k(s/)ds/ 


(25) 
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and 


Ms)  = 


- u,s0 

u(s)kg(s) 


,ke(sO,Pe(sO) 


r s _ 

/ c(s/)p(s/)k(s/)p  (s')ds' 

J 8.. 


(26) 


For  the  computation  of  the  radiance  value  that  includes  atmospheric  absorption, 
u(s0),  £e(s0),  and  (3e(s0)  are  computed  from  Eqs.  (10),  (7),  and  (8),  respectively, 
for  the  whole  atmospheric  path  by  using  s = Sq  as  the  upper  integration  limit. 

For  the  computation  of  the  unattenuated  source  radiance,  u(sQ)  is  simply  set  to 
zero.  This  corresponds  to  using  c(s')  = 0 in  Eq.  (10)  and  represents  a 
nonabsorbing  atmospheric  path. 

B.  EXACT  FORMULATION  FOR  VOIGT  LINE  SHAPE 

Formal  solutions  for  the  equivalent  width  derivative  functions  can  be 

obtained  for  the  Voigt  line  shape  by  the  same  procedures  used  earlier  for  the 

3 4 

Lorentz  and  Doppler  lines.  ’ One  has  only  to  carry  through  the  development 
with  the  spectral  absorption  coefficient  for  the  Voigt  line 

k(v,s)-^  -S^s)K[a(s),  w(s)]  (27) 

where  S(s)  is  the  lire  strength,  a(s)  is  a measure  of  the  local  Lorentz  line 
width  to  Doppler  line  width 


a(s)  = yin  2 


Yl  (s) 


(28) 
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w(s)  is  the  measure  of  spectral  displacement  from  the  line  center 


v — v r 


W(S)  = 


and  K(a,  w)  is  the  line  profile  function 


K(a,  w) 


2 

/■»  -u 

e 

a2  + (w  - 

00 


du 


U ) 


The  result  for  an  isolated  Voigt  line  is 


d = c(s)p(s)S(s)y[x(s),  a(s),  pL(s),  PD(s)] 


where 


c ® 

y(x,a,pL,  pD)  = ^ / K(a,  z)e 

Jr\ 


-xK(apD/pL>  zpD) 


dz 


The  arguments  of  y are 


x(s)  = 


/In  2 e 


S (s)u(s) 


" Ws) 


Yl(  = ) 


pL(a)  = veL(s| 
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-:a,. 


(20) 


(30) 


(31) 


(32) 


(33) 


(34) 


and 


yD(s) 

P°,S,=WS> 


(35) 


with  a(s)  given  by  Eq.  (2.8).  Note  that  the  Doppler  dimensionless  optical  depth 
is  used.  The  path-averaged  parameters  Se(s),  Ve^(s)>  an^  ^^(s)  were  defined 
earlier.  For  a band  of  equal  intensity  Voigt  lines 


^ ^ c(s)p(s)k(s)y[x(s),  a(s),  pT  (s),  p^s)] 


D' 


(36) 


where  y is  the  same  function  as  for  an  isolated  line,  Eq.  (32),  with  the  x and 
p arguments  defined  by  Eqs.  (5)  and  (6),  respectively,  and  a(s)  by 


| ?>r(s) 

a(s)  w F^T> 


(37) 


For  a band  of  lines  with  an  exponential  or  exponential -tailed  inverse  line 
strength  distribution,  the  mean  equivalent  width  derivative  is  given  by  Eq.  (36), 
but  the  derivative  function  y is  given  by 


, v 2 r°  K(a,  z) 

y X’  a,  PD*  Pl}  ' sfW  / ri  , Kl  , .in 

J0  l1  + ^K(apD/pL,  zpD)J 


dz 


*38) 


with  n = 2 for  the  exponential  distribution  and  n = 1 for  the  exponential -tailed 
inverse  distribution. 

None  of  the  expressions  for  y given  by  Eqs.  (32)  and  (38)  is  amenable  to 
a closed-form  solution.  For  the  pure  Doppler  line  shape,  solutions  were 
obtained  by  compiling  tables  of  y with  respect  to  the  two  independent  variables 
x and  p from  which  y could  be  obtained  by  interpolation.  For  the  Voigt  line 
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shape,  however,  there  are  four  independent  variables,  and  the  tabular  approach 
is  not  practical.  The  solution  used  here  is  to  approximate  the  equivalent  width 
of  a Voigt  line  by  a combination  of  the  component  Lorentz  and  Doppler  equivalent 
widths . 


C.  APPROXIMATE  FORMULATION  FOR  VOIGT  LINE  SHAPE 


The  Voigt  profile  is  complex  enough  that,  even  for  homogeneous  paths, 
simple  solutions  for  the  equivalent  width  of  an  isolated  line  have  not  been 
formulated.  The  equivalent  width  is  either  computed  numerically  as  needed, 
interpolated  from  tables  of  equivalent  width  versus  a and  x,  or  computed  from 
approximate  expressions.  An  example  of  available  interpolation  tables  has 
been  given  by  Jansson  and  Korb.  Three  methods  of  approximating  the  equiva- 
lent width  have  been  discussed  by  Hansen  and  McKenzie,  ^ Yamada,  ~ and 

g 

Rogers  and  Williams,  respectively.  The  method  of  Hansen  and  McKenzie 
approximates  W(a,x)  directly  ir.  terms  of  a ana  x.  Yamada's  approach  takes 
W to  be  the  Lorentz  or  Doppler  equivalent  width  in  the  limits  of  strong  and 
weak  absorption,  respectively,  and  gives  the  solution  in  the  transition  region 
as  one  or  the  other  of  these  equivalent  widths  modified  by  a function  of  a and 
x.  Of  the  three  approximations,  that  of  Rodgers  and  Williams  is  the  simplest 
to  implement  and  also  the  most  accurate.  For  a homogeneous  path  and  isolated 
Voigt  line,  their  representation  for  W is 


P.  A.  Jansson  and  C.  L.  Korb,  "A  Table  of  the  Equivalent  Width  of  Isolated 
Lines  with  Combined  Doppler  and  Collision  Broadened  Profiles,  " J.  Quant. 
Spectr.  Radiative  Transfer  8.  1399(1968). 


C.  F . Hansen  and  R.  L.  McKenzie,  "Analytic  Approximations  for  Curves  of 
Growth  of  Doppler -Lorentz  Broadened  Lines,"  J.  Quant.  Spectr.  Radiative 
Transfer  11,  349  (1971).  
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H,  Y.  Yamada,  "Total  Radiances  and  Equivalent  Widths  of  Isolated  Lines 
with  Combined  Doppler  and  Collision  Broadened  Profiles,  " J.  Quant.  Spectr. 
Radiative  Transfer  8.  1463  (1968). 


C.  D.  Rodgers  and  A.  P.  Williams,  "Integrated  Absorption  of  a Spectral  Line 
with  the  Voigt  Profile,  " J.  Quant.  Spectr.  Radiative  Transfer  14,  319  (1974). 
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r.  . 


2 2 

W + 

L D 


WLWD 


^2jl/2 


where  Wj  and  are  the  equivalent  widths  that  would  result  if  the  line  shape 
were  assumed  to  be  purely  Lorentz  and  Doppler,  respectively,  and  W^.  is  the 
weak-limit  form  for  the  equivalent  width  of  any  line  shape 


Ww  = Su 


The  reported  maximum  error  of  the  approximation  is  ~8  percent  and  occurs 
for  x (the  Doppler  depth)  when  a ^0.47  and  along  the  curve  ax  ^15  when 
a < 0.00083  (x  > 1.8  X 104). 

This  approximation  is  used  here  for  bands  of  lines  and  inhomogeneous 
optical  paths  as  follows.  The  mean  equivalent  widths  W^(s)/6and  W^(s)/6are 
computed  in  either  the  CG  or  LS  approximations  according  to  Eqs . (21)  or  (22) 
and  are  combined  according  to  the  algebraic  form  of  Eq.  (39)  in  order  to  obtain 
W(s)/6.  The  weak-limit  form  is  taken  as 


Ww(s) 


= u(s)ke(s) 


The  result  for  W(s)/6is  then  used  in  Eq.  (2)  to  obtain  the  transmittance  t(s). 
The  derivative  of  the  mean  equivalent  width  could  be  obtained  by  a numerical 
differentiation  of  W(s)/6,  but  is  obtained  in  practice  by  differentiating  Eq.  (39) 
analytically  to  obtain  the  derivative  function 


y(x,  a,  PD,  PL)  = AyL(x/2a  </tt,  PL)  + ByD(x,  pD)  + C 
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where  and  y arc  the  Lorcntz  and  Doppler  derivative  functions,  respectively. 
The  coefficients  of  Eq.  (42)  are 


WT  /6 

A - 

W/6 


M 

yv6/ 


(43) 


B = 


W^(6 

W/6 


\ww/s/ 


c 


ww/6  /wl/6\2/w-/6'2 


W/6 


D' 


W„,/6/\  W,.,/6 


W 


W' 


(44) 


(45) 


The  dimensionless  optical  depth  x is  taken  as  the  Doppler  depth.  With  this 
result  for  y,  the  mean  equivalent  width  derivative  is  computed  from  Eq.  (4) 
and  the  transmittance  derivative  from  Eq.  (3). 

The  main  assumptions  made  in  this  use  of  the  Rodgers  and  Williams 
approximation  is  that  the  mean  equivalent  widths  for  arrays  of  Lorentz  and 
Doppler  lines  combine  in  the  same  way  as  the  equivalent  widths  of  isolated 
lines.  Since  the  curve-of-growth  functions  for  bands  of  Lorentz  lines  do  not 
differ  greatly  from  that  for  an  isolated  line  regardless  of  the  line  intensity 
distribution  used  (and  similarly  for  the  Doppler  case),  this  assumption  should 
not  be  too  inaccurate. 

Note  that  the  functional  form  of  Eq.  (39)  always  yields  W /6  > W^/6, 

W/6  — Wj-^/6.  Consequently,  the  transmittance  t computed  by  Eq.  (2)  for  the 
Voigt  line  shape  is  always  less  than  or  equal  to  the  transmittance  computed 
for  the  pure  Lorentz  or  Doppler  line  shape.  This  result  is  not  necessarily 

true  for  the  effective  transmittance  t . 

e 
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III.  RESULTS  WITH  THE  CG  AND  LS  INHOMOGENEITY 

APPROXIMATIONS 


A.  CALCULATION  CONDITIONS  AND  DATA 

The  band  model  formulation  of  Section  II  is  complete,  but  its  use 
requires  specification  of  the  pressure,  temperature,  and  species  concen- 
tration (pTc)  variations  along  the  atmospheric  and  source  portions  of  the 
optical  path  and  specification  of  the  fundamental  band  model  parameters 

k,  6 , and  v. 
e 1 

The  pTc  variation  along  the  atmospheric  path  was  interpolated  from 
the  altitude  profiles  tabulated  by  McClatchey  et  al.^  for  five  model  atmo- 
spheres. For  all  of  the  results  presented  here,  the  tropical  model  atmo- 
sphere was  used.  The  profiles  for  this  atmosphere  are  plotted  in  Fig.  1 
The  line-of-sight  viewing  geometry  consists  of  a space  sensor  looking  down 
through  the  atmosphere  at  a hot  gas  plume  emission  source  at  some  alti- 
tude between  0 and  60  km.  The  zenith  angle  from  the  target  to  sensor  is 
75  deg. 

The  plume  model  used  is  given  in  Table  3.  The  model  assumes  a 
homogeneous  and  isothermal  hot  gas  exhaust  contained  in  an  axisymmetric 
cylindrical  volume.  The  pTc  properties  and  dimensions  of  the  cylinder  are 
allowed  to  vary  with  altitude  in  order  to  simulate  the  plume  properties  of  a 
real  missile.  All  of  the  results  obtained  here  pertain  to  a single  line  of 
sight  that  intersects  this  plume  at  right  angles  to  the  plume  axis  and  passes 
through  a full  diameter  of  the  plume.  The  source  portion  of  the  total  path 
is  thus  a homogeneous  and  isothermal  path  segment. 


9 

R.  A.  McClatchey,  R.  W.  Fenn,  J.  E.  A.  Selby,  F.  E.  Volz,  and 

J.  S.  Garing,  Optical  Properties  of  the  Atmosphere  (Revised),  AFCRL-71- 

0279,  Air  Force  Cambridge  Research  Laboratory,  Bedford,  Mass. 

(10  May  1971). 
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Preceding  page  blank 


Table  3.  Homogeneous  and  Isothermal  Plume  Model 


z,  km 

r,  m 

L,  m 

T,  °K 

0 

1.488 

16.21 

2207.0 

5 

1.718 

37.64 

1867.4 

10 

2.  170 

71.35 

1625. 2 

15 

2.  870 

133.9 

1417. 8 

20 

3.865 

237.9 

1249.4 

25 

5.  197 

409.  5 

1105.  1 

30 

6.897 

704.  1 

975.8 

35 

8.918 

1.260  X 103 

844.4 

40 

11.29 

2.237  x 103 

723.5 

45 

14.09 

3.789  X 103 

622.5 

50 

17.54 

5.944  x 103 

545.6 

55 

22.46 

9. 817  x 103 

478.7 

60 

29.41 

1.444  x 104 

441 . 5 

z = source  altitude;  r = plume  radius;  L = plume  length;  T = 
plume  temperature.  The  plume  gas  pressure  is  the  ambient 
pressure  at  altitude  z.  The  mole  fraction  of  H2O  is  0.  2533 
for  all  z.  The  mole  fraction  of  C02  is  0.03714  for  all  z. 


The  band  model  parameters  used  here  were  discussed  in  detail  in  an 
earlier  report.  For  both  and  I^O,  the  parameters  represent  a com- 
bining of  band  model  parameters  derived  from  the  comprehensive  atmospheric 


S.  J.  Young,  Band  Model  Parameters  for  the  2.7-um  Bands  of  H?Q  and 
CO?  in  the  100  to  300°K  Temperature  Range.  TR-0076(6970) -4,  The 
Aerospace  Corporation,  El  Segundo,  Calif.  (31  July  1975). 
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absorption  line  data  compilation  of  AFCRL**  and  band  model  parameters 

tabulated  in  the  NASA  Handbook  of  Infrared  Radiation  from  Combustion 
12 

Gases.  The  former  band  model  parameters  were  used  to  represent  the 
low-temperature  variation  of  the  final  parameter  set;  the  latter  represent 
the  high-temperature  variation.  For  intermediate  temperatures,  the  param- 
eters of  the  final  set  were  interpolated  from  transition  curves  drawn  from 
one  of  the  component  sets  to  the  other.  For  both  C02  and  H2Ot  the  final 
parameter  sets  display  a smooth  and  consistent  variation  of  the  parameters 
from  temperatures  as  low  as  100  to  as  high  as  3000  °K.  The  sets  cover  the 
entire  2.  7 -pm  bands  of  H2<D  and  C02  and  were  constructed  to  represent  a 
spectral  resolution  of  25  cm”*. 

The  line  spacing  parameter  of  the  final  set  is  consistent  for  use  with 
the  statistical  band  model  for  a constant  line  intensity  distribution.  For  other 
distributions,  6^  must  be  multiplied  by  an  appropriate  constant  or  this  con- 
stant must  be  incorporated  in  the  equations  of  the  band  model  formulation. 

The  latter  approach  is  taken  here  and  explains  the  factors  of  tt  (the  appro- 
priate constant  for  the  exponential-tailed  inverse  line  strength  distribution) 
that  occur  at  various  places  throughout  the  discussion  of  Section  II  (and  later 
in  Section  IV). 

Most  of  the  results  presented  here  are  for  the  transmittance  7 and  the 
effective  transmittance  (each  computed  with  25  cm  * resolution)  appropri- 
ately averaged  over  three  wide  bandpasses  in  the  2. 7 -pm  region.  The  re- 
gions were  selected  to  represent  the  low  wavenumber  wing  (3350-3575  cm-*), 
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the  bandcenter  (3575-3725  cm’*),  and  the  high  wavenumber  wing  (3725-3950 
cm’*)  of  the  combined  W^O/CO^  vibration  band.  Radiative  transfer  effects 
in  the  latter  spectral  region  are  dominated  by  F^O.  If  v^  and  v ^ represent 
the  lower  and  upper  wavenumber  limits,  respectively,  of  one  of  these  band- 
passes,  the  average  transmittances  for  t and  t are  computed  from 


Jv . 


t(v)  dv 


L (v)  dv 
a 


T = — - 
e /.v 


•'M 


(v)  dv 


where  L (v)  and  L (v)  are  the  25  cm  resolution  solutions  of  Eq.  (23)  for  the 
a u 

cases  where  the  intervening  atmosphere  is  and  is  not  absorbing,  respectively 

Since  L (v)  = t (v)  L (v)  by  definition,  the  form  of  Eq.  (47)  reduces  to  that  of 
a e u 7 

Eq.  (46)  if  Lu(v)  is  constant  across  the  bandpass. 

B.  RESULTS  AND  DISCUSSIONS 


Bandpass  averaged  results  for  T versus  altitude  are  shown  in  Fig.  2. 
For  all  three  bandpasses,  the  curve  for  the  Voigt  line  shape  makes  a smooth 
transition  from  the  Lorentz  curve  to  the  Doppler  curve  in  the  altitude  region 
from  15  to  40  km.  The  spectral  variations  of  t (25  cm  * resolution)  for  the 
altitude  z = 35  km  and  for  each  line  shape  are  shown  in  Fig.  3,  where  it  is 
shown  that,  at  all  spectral  positions,  the  absorption  caused  by  the  Voigt  line 
shape  is  the  greatest. 


Figure  2.  Altitude  Variation  of  T Computed  for  the 
Lorentz,  Doppler,  and  Voigt  Line  Shapes 
and  with  the  CG  Approximation,  (a)  3350- 
3575  cm-*  bandpass;  (b)  3575-3725  cm'* 
bandpass;  (c)  3725-3950  cm-*  bandpass. 
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All  of  the  results  or  Figs.  2 and  3 were  obtained  by  using  the  CC.  ap- 
proximation. The  results  obtained  with  the  LS  approximation  were  nearly 
the  same  as  with  the  CG  approximation.  The  largest  difference  in  the  result 
for  T computed  with  the  two  approximations  was  - 60  percent,  but  this  oc- 
curred for  the  center  bandpass  (strongest  absorbing  region)  and  z , 0.  where 
t was  ~ 2.4  x 10  . The  largest  difference  for  a transmittance  value  large 

enough  to  be  meaningful  was  ~ 7 percent  and  occurred  for  the  center  bandpass 
with  z = 10  km.  The  transmittance  value  was-  0.08.  Above  - 15  km,  the 
difference  in  results  obtained  with  the  CG  and  LS  approximations  was  5 1 per- 
cent, and  the  transmittance  s were  all  > 0.3.  These  results  reinforce  the  hy- 
pothesis that  the  degree  of  inhomogeneity  encountered  along  a strictly  atmo- 
spheric slant  path  is  not  too  large  to  be  adequately  handled  by  the  CG 
approximation. 

Bandpass  averaged  results  for  Te  versus  target  source  altitude  are 
shown  in  Figs.  4 and  5 for  calculations  made  with  the  CG  and  LS  approxima- 
tions, respectively.  Comparison  of  results  is  facilitated  by  discussing  the 
results  for  each  line  shape  individually.  Consider  first  the  Lorontz  line  shape 
results.  The  failing  of  the  CG  approximation  for  highly  inhomogeneous  opti- 
cal paths  and  Lorentz  line  shape  has  been  shown  to  result  in  computed  values 
for  Te(v)  that  are  too  large  at  some  spectral  positions;1  so  large,  in  fact, 
that  frequently  t^v)  is  not  only  greater  than  t(v),  but  also  greater  than  unity. 
These  errors  occurred  in  regions  of  weak  absorption,  such  as  the  band  wings. 
An  example  of  this  failing  is  demonstrated  in  Fig.  6a  for  the  altitude  z = 15  km. 
The  LS  approximation  was  devised  to  correct  this  discrepancy,  and  the  result 
of  Fig.  6b  displays  its  success.  Although  tJv)  > 1 for  some  v in  the  CG  ap- 
proximation, the  extent  of  the  error  is  not  large  enough  to  make  itself  felt  in 
the  bandpass  averages  T^.  The  Lorentz  line  results  for  T in  Fig.  4 are  all 

less  than  unity.  They  are,  however,  greater  than  the  LS  results  of  Fig.  5 in 
every  case. 
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Figure  4.  Altitude  Variation  of  Te  Computed  for  the  Lorentz,  Doppler 
and  Voigt  Line  Shapes  and  with  the  CG  Approximation. 

(a)  3350-3575  cm"*  bandpass;  (b)  3575-3725  cm"*  bandpass; 
(c)  3725-3950  cm"*  bandpass. 
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Figure  5.  Altitude  Variation  of  Te  Computed  for  the  Lorentz, 
Doppler,  and  Voigt  Line  Shapes  and  with  the  LS 
Approximation,  (a)  3350-3575  cm-1  bandpass; 

(b)  357503725  cm-*  bandpass;  (c)  3725-3950  cm"* 
bandpass . 
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Figure  6.  Spectral  Variation  of  T (v)  at  z = 15  km 
Computed  for  the  Lorentz  Line  Shape. 

(a)  CG  approximation  [Te(v)  is  not  phys- 
ically meaningful  whenever  it  is  greater 
than  t(v)]j  (b)  LS  approximation. 
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In  almost  every  case  of  Figs.  4 and  5.  the  result  for  T obtained  for 
the  pure  Lorentz  line  shape  increases  with  z up  to  ~ 15  to  20  km.  and  then 
maintains  a value  much  less  than  unity  as  z is  increased  further.  This  ef- 
fect was  mentioned  in  the  Introduction  and  is  a result  of  using  the  pressure  - 
broadened  Lorentz  line  shape  in  an  altitude  region  where  its  use  is  not  valid 

The  Tc  results  for  the  Doppler  line  shape  computed  with  the  CG  ap- 
proximation are  shown  in  Fig.  4.  As  for  the  Lorentz  line  shape,  7 (v)  com- 
puted with  this  approximation  can  be  greater  than  unity.  Now,  however,  the 
seriousness  of  the  error  is  so  great  that  even  the  bandpass  averages  can  ex- 
ceed unity.  For  the  3350  to  3575  cm-1  bandpass  (Fig.  4a),  T begins  to  ap- 
proach a realistic  variation  with  altitude  at  z - 30  km.  For  the  o.her  two 
bandpasses,  a realistic  variation  is  not  yet  achieved,  even  at  z = 60  km.  For 
the  Lorentz  line  case,  the  LS  approximation  was  able  to  give  significantly 
more  realistic  results  for  Tjv)  than  did  the  CG  approximation.  The  results 
of  Fig.  5 indicate  that,  for  Doppler  lines,  the  LS  approximation  gives  only 
little  improvement  over  the  CG  approximation.  Values  of  T >1  occur  for 
each  of  the  bandpasses.  The  slight  improvement  provided  by  the  LS  approxi- 
mation is  that  the  variation  of  Te  with  altitude  becomes  reasonable  at  a lower 
altitude.  For  the  low  wavenumber  bandpass,  the  variation  is  reasonable 

above  z ~ 20  km.  For  the  other  two  bandpasses,  the  variation  is  reasonable 
above  z ~ 40  km. 

The  results  for  Te  obtained  for  the  Voigt  line  shape  are  also  shown  in 
Figs.  4 and  5.  Because  the  calculations  for  the  Voigt  line  shape  are  carried 
out  with  the  results  for  the  pure  Lorentz  and  Doppler  line  shapes,  and  be- 
cause the  Doppler  component  results  for  both  the  CG  and  LS  approximations 
are  unrealistic,  the  results  for  the  Voigt  line  shape  do  not  have  much  mean- 
ing. Some  interesting  features  should  be  mentioned  however.  The  Voigt  re- 
sults for  Te  computed  with  the  LS  approximation  always  lie  between  the  T 
curves  computed  for  the  Lorentz  and  Doppler  line  shape.  With  the  CG  ap- 
proximation, on  the  other  hand,  this  behavior  is  not  obtained.  In  fact,  the 
results  obtained  with  the  CG  approximation  display  a surprisingly  realistic 


mm 


• variation  for  computed  for  the  Voigt  line  shape  despite  the  unreasonableness 

of  the  component  line  shape  results.  The  behavior  of  neither  the  CG  nor  the 
LS  approximations  follows  that  required  for  T;  viz.,  that  the  Voigt  result  for 
T be  smaller  than  either  of  the  component  line  shape  results. 

| 

In  the  next  section,  the  failure  of  the  LS  approximation  when  used  with 
the  Doppler  line  shape  is  discussed,  and  a new  approximation  is  formulated  to 


account  for  inhomogeneities  along  an  optical  path. 


IV.  THE  TWO-PATH  DERIVATIVE  APPROXIMATION 


A.  INTRODUCTION 

In  order  to  understand  the  fundamental  reasons  for  the  failure  of  the  LS 
approximation  when  applied  to  the  Doppler  line  shape,  consider  the  simplified 
problem  of  two  homogeneous  path  segments  in  series.  The  first  lies  between 
s - 0 and  s = s^  and  has  optical  and  thermodynamic  properties  identified  by 
subscript  1;  the  second  lies  between  s = sQ  and  s = S and  has  properties  iden- 
tified by  subscript  2.  Two  effects  must  be  considered:  the  change  in  line 
width  as  s traverses  from  the  first  path  segment  into  the  second,  and  the 
change  in  line  density  for  the  same  transition.  Consider  first  an  isolated  line, 
and  in  order  to  simplify  the  problem  even  further,  assume  the  line  to  be  rec- 
tangular in  shape.  The  variation  of  the  line  width  with  s is 


Y(s)  = 


Yj  0<  s < sQ 

y2  sq  < s < S 


(48) 


Assume  that  y2  > Yj  • (As  an  approximation  for  Doppler  lines,  y2  > Yj  im- 
plies that  the  temperature  relation  for  the  two  paths  is  T2  > Tj.  For  JLorentz 
lines,  Y2  > Yj  would  occur,  for  example,  if  the  pressure  relation  were 
P2  > Pj  when  all  the  other  properties  of  the  two  path  segments  were  equal.  ) 
The  variation  of  the  path -averaged  line  width  with  s is 


0 < s < s, 


Ye(s)  = { 


:lPl^lSlsQ  + C2P2V2S2  (S  ~ V 
ClplSls0  + C2P2S2  (s  - S0)  sQ<s<S 


(49) 
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The  variations  of  y(s)  ar>d  Ye(s)  3-re  shown  schematically  in  Fig.  7 for  a 
rectangular  line  centered  on  Vq.  In  the  LS  approximation,  the  contribution 
to  the  radiance  at  s = 0 from  a path  element  ds  at  s is  obtained  by  assuming 
that  the  emission  line  width  in  ds  is  y(s),  and  that  the  radiation  from  ds  is 
attenuated  by  a homogeneous  path  between  s = 0 and  s along  which  the  absorp- 
tion line  has  the  constant  value  Ye(s)-  As  is  evident  from  Fig.  7,  when 
S ">  S0’  ^e(s)  < Y(8>.  Consequently,  the  radiation  emitted  in  ds  from 
any  spectral  element  dv  in  the  line  wings  beyond  | v - Vq  | - Y (s)  is  not 
attentuated  at  all  in  traversing  the  path  to  s = 0 . If  the  path  properties  are 
such  that  y2  » v i’  and  the  path-averaged  line  width  Ye(s)  is  so  heavily 
weighted  by  the  properties  of  region  1 that  Ye(s)  — Yj>  even  for  s >>  Sq,  a 
substantial  part  of  the  radiation  from  region  2 will  arrive  at  s = 0 without 
being  attenuated.  On  the  other  hand,  when  the  radiance  at  s = 0 due  to  the 
second  region  alone  is  computed,  no  consideration  is  given  to  region  1,  and 
Ye( s ) has  the  constant  value  Ye(s)  = V2  -or  a11  s > sQ.  In  this  case,  the  radia- 
tion from  each  spectral  interval  dv  of  ds  is  properly  attenuated  as  it  traverses 
toward  s = 0.  The  net  effect  is  that  a larger  radiance  is  computed  for  the 
total  path,  which  consists  of  regions  1 and  2 in  series,  than  is  computed  for 
region  2 alone.  This  result  is  totally  unacceptable . For  real  spectral  line 
shapes,  the  wing  radiation  from  the  emission  line  is  not  completely  unattenu- 
ated in  the  total  path  radiance  calculation,  but  can  be  significantly  underesti- 
mated because  of  this  effect. 

The  consideration  of  the  change  in  line  density  between  the  two  regions 
is  similar  to  that  for  the  line  width.  Assume  that  Yj  = Y2  anc*  T2  > Tj.  Then, 
the  line  density  (1/6  ) will  be  larger  in  region  2 than  in  region  1.  and  3,  >3,. 
The  path  variations  of  3(s)  and  3e(s)  will  be  given  in  Eqs.  (48)  and  (49),  but 
with  y's  replaced  by  3's-  Again,  in  the  LS  approximation,  the  emission  from 
ds  is  assumed  to  be  described  by  the  line  density  parameter  3(s)  and  the  at- 
tenuation by  a homogeneous  path  with  line  density  parameter  3p(s)-  Clearly, 
3e(s)  < 3(s)  for  s > Sq.  In  this  case,  the  radiation  from  complete  lines 
(rather  than  from  just  line  wings)  will  be  improperly  attenuated  in  the  total 
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path  radiance  computation  because  each  line  in  region  2 will  not  have  a 
corresponding  line  in  the  attenuation  path  to  s — 0 that  can  absorb  its  emission 
radiation.  For  region  2 alone,  P^(s)  = p^  and  each  emission  line  does  have  a 
corresponding  absorption  line. 

The  limitation  of  the  LS  approximation  is  that,  for  optical  paths  along 
which  y(s)  or  P(s)  increases  sharply  at  some  point  (sq).  the  path-averaged 
quantities  y^(s)  or  3^(s)  do  not  respond  to  the  change  fast  enough  to  give  values 
sufficiently  large  to  properly  account  for  the  attenuation  of  radiation  from  path 
elements  beyond  the  region  of  sharp  increase.  The  solution  to  this  limitation 
is  to  allow  the  attenuation  of  radiation  from  ds  for  s > sQ  to  be  effected  by  two 
homogeneous  paths  that  lie  between  s = 0 and  s.  The  first  path  is  the  whole  of 
region  1 (s  - 0 to  s = Sq),  and  the  absorption  line  width  or  density  parameter 
is  taken  as  the  average  value  for  this  whole  region.  The  second  path  is  the 
portion  of  region  2 between  s = sQ  and  s,  and  the  line  width  or  density  param- 
eter is  the  average  value  between  sQ  and  s.  The  effect  is  that  yp(s)  or  p (s) 
is  able  to  follow  the  variations  of  y(s)  or  p(s)  at  s - sQ  because  their  compu- 
tation is  begun  anew  at  s s^  and  is  not  influenced  by  the  properties  of  the 
path  for  s < sQ.  The  formulation  of  this  approach  for  the  Lorentz  and  Dop- 
pler line  shape  is  given  in  the  following  sections. 

Because  the  formulations  involve  two  homogeneous  paths  and  derive 
expressions  for  the  derivative  of  the  equivalent  width,  the  approximation  is 
referred  to  here  as  the  two-path  derivative  approximation.  The  LS  approxi- 
mation corresponds  to  a single -path  derivative  approximation.  Similarly, 
the  present  use  of  the  CG  approximation  could  be  considered  as  a single -path 
equivalent  width  approximation.  The  corresponding  two-path  equivalent  width 
approximation  would  be  the  models  presented  by  Plass.  14 


G.  N.  Plass,  "Spectral  Band  Absorptance  of  Radiation  Traversing  Two  or 
More  Cells  in  Series,  " Appl,  Opt.  4,  69  (1965). 

G.  N.  Plass,  "The  Theory  of  the  Absorption  of  Flame  Radiation  by  Molecular 
Bands,  " Appl.  Opt.  4,  1 6 1 (1965). 


B.  FORMULATION  FOR  LORENTZ  LINE  SHAPE 

The  equivalent  width  derivative  for  an  isolated  Lorentz  line  is 


ds 


r 

fB 

I k(v,  s)  exp 

- I c(s'  )p(s')k(v,  s')  ds' 

/-CO 

Jo 

dv  (50) 


where  the  spectral  absorption  coefficient  is 


Mv.  <0  = 3 1} XlL 

(v  - v„r  + v (») 


(51) 


The  exponential  factor  of  Eq.  (50)  describes  the  attenuation  of  radiation  from 
the  path  element  ds  by  the  path  between  s = 0 and  s,  and  it  is  in  the  path  inte- 


gral of  this  factor  that  the  two-path  approximation  is  made.  For  s > sQ,  the 


path  integral  can  be  written  [by  using  Eq.  (51)]  as 


f c(B')p(s')k(v,  s')ds/  = - f 0 £.(8/)P(8/)S(8/)y(s/)  ds/ 
JO  ^ J0  ( v-v  ) +/(s') 


+ 1 f c(S')p(S')S(s')y(s'] 
ir  I . .2  2 

Js  (V  " v0)  + Y (s/) 


(52) 


For  the  integral  between  s = 0 and  s = Sq,  y(s7)  in  the  denominator  of  the  inte- 


grand is  replaced  by  the  path-averaged  value  between  s = 0 and  s = sQ, 


Vl(s0)-  Similarly,  for  the  second  integral,  y(s')  is  replaced  by  the  average 
between  s = sQ  and  s,  V2(8)-  Then, 
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Substitution  of  Eq.  (53)  into  (50),  the  change  of  variable  z = 
the  use  of  the  definitions 


(v  " vQ)/ y(s),  and 


p.(8)  = 

1 v,(V 


(62) 


and 


(63) 


yield 


dW(s)  _ 

d8  - c(s)  p(s)  S(s)  y[  Xl(s0),  x2(s),  Pl(s),  p2(s)] 


(64) 


where  the  equivalent  width  derivative  function  y(x.,  x_,  p , p ) 

12  12 


is 


y(Xl’  x2’  pr  p2*  = n 


lvo  1 

2xl  &2 
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As  in  the  LS  approximation,  a closed-form  solution  for  Eq.  ,65)  could  no,  to 
lound.  However,  for  bands  of  lines  in  which  the  line  strengths  are  distributed 
at  < ording  to  an  exponential  or  exponential -tailed  inverse  line  strength  distri- 
button,  closed-form  solution,  can  be  obtained  by  tntegrat.ng  in  the  complex  a 
plane  by  ustng  the  method  of  residues.  Arguments  presented  earlier  3 ' 4 allow 
the  mean  equivalent  width  derivative  for  bands  of  lines  to  be  derived  as 

1 dW  ( s ) 

6 ds  cts)  p(s)k(sl  y[oc,(s0),  x2(s),  Pjls),  p ,,)]  (66) 
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and  wheie  n 2 for  the  exponential  line  strength  distribution  and  n 1 for  the 
exponential-tailed  inverse  distribution.  The  factors  that  account  for  com. 
patabi lity  with  the  band  mode,  parameter  6,  are  suppressed  in  Eq . ,67.  but 
must  be  included  in  the  final  solutions  for  y.  For  n 2.  both  x,  and  x on 
the  right-hand  stde  of  Eq.  (67)  must  be  multiplied  by  s/4:  for  n I.  they  must 

be  multiplied  by  it.  For  bands  of  lines,  the  x and  p parameters  of  Eq.  [661 
are  defined  by 


ki'V  W 


W = / c(s/)  p(s/)  ds/ 


-44- 


(S0>  =u"1(8q)  J c(8')  p(s')k(6')ds' 


. r° 

Z / C(s', 

ui<8o>ki,so,-'o 


) p(s')  k(s')  g(s')  ds 


k2(s)  u2(s) 


x2,8>  = 32(.)- 


U2(S)  =j  C(8')p(s') 


ds 


(8) 


OT)  f c(e<)  P(S')  k(8' 


) ds' 


i /*8 

Z / c(s')  p(s')  k(s')  0(s')  ds' 

u2(8)k2(s)7s 


pj(s)  = ^(sJ/P^Sq) 


p9(s)  -•  P(s)/p,(s) 


A straightforward,  but  tedious.  manipulation  of  Eq . (67)  allows  it 


be  written  as 
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The  integration  of  Eq.  (78)  by  the  method  of  residues  is  accomplished 

by  closing  the  contour  from  to  a by  an  infinite  semicircle  in  the  upper  half 

of  the  complex  z plane.  Examination  of  the  integrand  of  Eq.  (7  8)  reveals  that 

the  only  poles  enclosed  by  the  contour  are  those  at  z^,  z y and  z^.  The  pole 

at  z.j  is  first  order,  whereas  those  at  z^  and  z^  are  iirst  order  for  n = 1 and 

second  order  for  n - 2.  The  calculation  of  the  residues  at  z.,  z~,  and  z,  bv 
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where  I(z)  denotes  the  integrand  of  Eq.  (78),  and  m is  the  order  of  the  pole  at 


z y is,  again,  straightforward  but  tedious,  and  only  the  final  result  is  given. 
For  n = 1 or  2 
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(5) 

and  a ^ is  of  the  same  form,  but  with  p^  and  p^  interchanged.  For  n = 2 
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(5) 

and  a _j  is  obtained  by  interchanging  Pj  and  p2>  The  final  formal  solution  for 
y(x  j , x2>  pJf  P2)  is,  then 


/ . 2 tt  i 

y*xl'  x2'  pl*  p2  " 2 2 

»P,  P2 

the  2i  factor  of  Eq.  (88)  cancels  with  the  2i  factors  in  the  denominators  of  the 
a_t  terms.  It  can  be  verified  that,  for  Xj  = 0,  n = i,  and  x2  = ttx,  this  solu- 
tion reduces  to  the  LS  approximation  solution  of  Eq . (13). 

C.  FORMULATION  FOR  DOPPLER  LINE  SHAPE 

The  derivation  of  the  y(xJ(  x2>  p p2)  functions  for  a Doppler  line  or 
bands  of  Doppler  lines  follows  that  for  the  Lorentz  line,  except  that  the  Dop- 
pler line  shape 
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-1 
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is  used.  The  formal  solution  for  an  isolated  line  is 


i 


t 


i 

l 


I 


For  the  exponent. al  (n  - 2)  and  exponential-tailed  inverse  (n  = 1)  line  strength 
distributions 


The  definitions  for  the  path  averages  are  the  same  as  for  the  Lorentz  case. 

Eqs  . (54)  through  (63)  and  (68)  through  (77),  except  that  Xj  and  x^  are  given  by 
equations  of  the  form 

J 


for  the  insolated  line  case,  Eqs.  (54)  and  (58),  respectively.  The  factors 
rr/4  for  n = 2 and  rr  for  n = 1 have  been  suppressed  in  Eq.  (91). 

Neither  Eq.  (90)  nor  (91)  is  amenable  to  closed-form  solutions.  In  the 
LS  approximation,  tables  of  values  were  prepared  from  which  the  desired 
value  could  be  interpolated.  However,  even  this  method  is  not  practical  here 
because  there  are  now  four  independent  variables  (Xj,  x^.  p j,  and  p^)  rather 
than  just  the  two  (x,  p)  of  the  LS  approximation  on  which  to  interpolate.  Tables 
for  adequate  ranges  of  the  four  variables  would  be  prohibitively  voluminous. 

The  solution  used  here  is  to  approximate  the  Doppler  line  shape  by  a parabolic 
line  shape.  This  approximation  is  applied  and  discussed  in  the  following 
section. 
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D. 


FORMULATION  FOR  PARABOLIC  LINE  SHAPE 


The  spectral  absorption  coefficient  for  a parabolic  line  shape  is 


k(  v,  s) 
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4^/2  V(  = ) 
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v ~ vq  \<  /Z  y(s) 
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(92) 


In  order  to  approximate  the  Doppler  line  shape  by  the  parabolic  line  shape, 
the  line  strength  S(s)  and  the  value  of  k(v,  s)  at  the  line  center  are  chosen  to 
be  the  same  for  the  two  profiles.  This  requires  that  the  parabolic  line  width 
be  chosen  as 


= 1. 127  y 


D 


(93) 


where  Hoppler  line  width.  The  resulting  profile  approximation  is 

shown  in  Fig.  8 

The  rationale  for  using  the  parabolic  approximation  is  threefold.  First, 
for  many  applications,  account  of  Doppler  broadening  provides  only  a small 
correction  to  the  radiative  transfer  effects  computed  with  the  Lorentz  line 
shape.  The  error  inherent  in  using  a parabolic  rather  than  Doppler  profile 
to  compute  this  small  correction  term  is  thus  minimized.  Second,  for  many 
cases  where  the  Doppler  effect  is  the  dominant  line -broadening  mechanism, 
the  absorption  for  the  path  is  small.  For  weak  absorption,  radiative  transfer 
effects  are  independent  of  line  shape  and  depend  only  on  the  line  strength. 
Third,  closed-form  solutions  can  be  obtained  for  the  parabolic  line  shape. 

Following  the  procedure  for  the  Lorentz  line,  the  equivalent  width 
derivative  function  for  an  isolated  parabolic  line  is 
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dz 


for  isolated  lines.  In  all  cases,  y is  the  parabolic  line  width.  Again,  the 
tt/4  and  it  factors  for  n = 2 and  1,  respectively,  are  suppressed  in  Eq.  (95). 

The  infinite  integration  limits  of  Eqs.  (94)  and  (95)  are  formal  limits 
only.  Because  the  parabolic  line  shape  has  a discontinuous  slope  at  the  point 
where  it  goes  to  zero,  the  integration  of  Eqs.  (94)  or  (95)  must  be  handled 
as  a set  of  six  cases,  which  result  from  the  six  inequality  relations  that  can 
be  formed  among  the  three  line  widths  Yj(sq).  Y2(s).  and  Y(s)-  Other  than 
the  tedium  of  keeping  track  of  these  cases,  the  integrations  are  all  elemen- 
tary, and  the  result  for  all  cases  can  be  written  as 
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For  an  isolated  line,  or  band  of  equal  intensity  lines 
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where  x xjPj  + x^p^,  an<^  Dawson's  function 
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For  a band  of  parabolic  lines  with  an  exponential  line  strength  distribut 
(n  = 2) 
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where  o'  = 1 + x j + x^. 

For  the  exponential-tailed  inverse  distribution 
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E.  RESULTS  AND  DISCUSSION 

The  calculations  with  the  two-path  derivative  approximation  were 
performed  exactly  as  for  the  CG  and  LS  approximations,  except  that  for 
s > sQ,  y was  computed  according  to  Eq.  (88)  rather  than  Eqs.  (11)  or 
(13)  for  the  Lorentz  component  of  the  line  and  to  Eqs.  (96)  and  (102)  rather 
than  Eq.  (16)  or  Table  1 for  the  Doppler  (parabolic)  component.  Also,  for 
S ^ s0’  tke  weak-limit  form  of  the  equivalent  width,  required  in  the  com- 
bining of  the  Lorentz  and  Doppler  (parabolic)  components  to  give  the  Voigt 
line  result,  was  taken  as 


Ww(s) 


VV  ki*so)  + U2(S)  k2(s) 


( 104) 


The  variation  of  Tg  with  altitude  for  the  three  bandpasses  is  shown  in 
Fig.  9.  The  results  for  the  Lorentz  line  component  are  indistinguishable 
from  those  obtained  with  the  LS  approximation.  The  largest  difference  that 
occurred  between  T^  computed  with  the  two  approximations  was  ~ 4 percent. 
The  results  for  the  Doppler  line  component  show  a substantial  improvement 
over  those  obtained  with  either  the  CG  or  LS  approximations.  The  variation 
with  altitude  is  physically  realistic  for  the  whole  altitude  region.  Spectra 
Te(v)  for  the  Doppler  line  shape  and  z = 30  km  are  shown  in  Fig.  10  for  all 
three  inhomogeneity  approximations. 

For  all  three  bandpasses,  the  result  for  the  Voigt  line  shape  displays 
a smooth  transition  from  the  Lorentz  curve  at  z 2;  15  km  to  the  Doppler  curve 
at  z ~ 50  km,  although  only  for  the  low  wavenumber  bandpass  is  the  transition 
to  the  Doppler  curve  complete.  As  for  the  LS  approximation,  the  Voigt  re- 
sult always  lies  between  the  Lorentz  and  Doppler  results. 


Spectra  of  Tg(v)  computed  for  each  line  shape  and  for  z = 35  km  are 
presented  in  Fig.  11.  To  within  calculational  accuracy,  the  Voigt  line  shape 
spectrum  always  lies  between  the  spectra  for  the  Lorentz  and  Doppler  line 
components  or,  in  the  band  wing  regions,  is  the  same  as  the  Doppler  spectrum. 
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F igure  9.  Altitude  Variation  of  Te  Computed  for  the  Lorentz, 

Doppler,  and  Voigt  Line  Shapes  and  with  the  Two-Path 
Derivative  Approximation,  (a)  3350-3575  cm'*  band- 
pass; (b)  3575-3725  cm'*  bandpass;  (c)  3725-3950  cm"* 
bandpass . 
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Figure  10.  Spectral  Variation  of  t(v)  at  z = 30  km  Computed  for  the 
Doppler  Line  Shape  and  with  the  CG,  LS,  and  Two-Path 
Derivative  Approximations.  For  the  latter  approxima- 
tion, the  line  shape  is  assumed  to  be  parabolic. 
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Figure  11.  Spectral  Variation  of  Te(v)  at  z = 35  km  Computed 
for  the  Lorentz,  Doppler,  and  Voigt  Line  Shapes 
and  with  the  Two-Path  Derivative  Approximation 


r 


,.  The  two-path  derivative  approximation  provides  a significant 

improvement  over  the  LS  approximation  for  the  Doppler  line  shape,  but  not 
for  the  Lorentz  line  shape,  at  least  for  the  problem  considered  here.  This 
effect  is  most  probably  a result  of  the  relative  shapes  of  the  two  line  pro- 
files. The  Doppler  line  has  little  wing  structure  and  more  closely  resembles 
a rectangular  line  (for  which  the  LS  approximation  is  worst)  than  does  the 
Lorentz  line  shape. 

In  Fig.  12,  the  3575-3725  cm"  1 bandpass  result  for  Tp  computed  for 
the  Voigt  line  shape  with  the  two-path  derivative  approximation  is  shown  to- 
gether with  the  atmospheric  transmittance  T for  the  same  line  shape  and 
bandpass.  The  transmittance  T is  the  correct  transmittance  to  use  for  com- 
putation of  the  attenuation  that  continuum  radiation  (such  as  black  body  radia- 
tion) or  line  emission  from  any  molecular  species  other  than  CC>2  or  H^O 
' „ would  experience  in  traversing  the  atmospheric  path  from  the  altitude  z to 

space.  For  continuum  radiation,  line  correlation  is  a meaningless  concept, 
whereas  for  non-H^O  or  CO2  radiation,  the  line  arrangement  between  the 
emission  spectrum  and  the  absorption  spectrum  is  almost  certain  to  be  un- 
correlated. This  transmittance  approaches  unity  rapidly  and  is  greater  than 
0.95  for  all  z > 40  km.  For  all  continuum  radiation  sources  (or  non-H?0/C0_ 
gaseous  emitters)  above  ~ 40  km,  the  atmosphere  can  thus  be  considered  as 
transparent.  For  the  hot  I-^O /CO^  plume  model  source  used  here,  the  cor- 
rect attenuation  factor  is  T . Its  value  of  ~ 0.7  at  40  km  can  hardly  be  con- 
sidered as  representative  of  a transparent  atmosphere.  Even  at  the  maximum 
altitude  considered  here  (z  = 60  km),  Tg  has  obtained  a value  of  only  ~ 0.  83. 
At  60  km,  on  the  other  hand,  T is  nearly  0.99.  Now  that  Doppler  broadening 
effects  I ave  been  realistically  accounted  for  by  using  the  two-path  derivative 
approximation,  this  remaining  large  difference  between  T and  T£  can  be  at- 
tributed only  to  line  correlation  effects  that  exist  between  the  hot  HgO/COg 
emission  spectrum  and  he  cool  H2O/CO2  atmospheric  absorption  spectrum. 
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Figure  12.  Altitude  Variation  of  T and  Te  for  the  3575- 
37  25  cm'*  Bandpass  Computed  for  the  Voigt 
Line  Shape  and  with  the  Two-Path  Derivative 
k Approximation 


LABORATORY  OPERATIONS 

Th«  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting 
experimental  and  theoretical  investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver- 
satility and  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel- 
opments is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratc  Ties  that  contribute  to  this  research  are: 

Aerophysica  Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans- 
fer, reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high-power  gas  lasers. 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  atmos- 
pheric optics,  chemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
1***r"i“d“ced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo- 
sensitive  materials  and  sensors,  high  precision  laser  ranging,  and  the  appli- 
cation of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics;  quantum  electronics, 
lasers,  and  electro-optica;  communication  sciences,  applied  electronics,  semi- 
conducting, superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory:  Development  of  new  materials;  metal 
matrix  composites  and  new  forms  of  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment;  application  of  fracture  mechanics  to  stress  cor- 
rosion and  fatigue-induced  fractures  in  structural  metals. 

Space  Physics  Laboratory:  Atmospheric  and  ionospheric  physics,  radia- 
tion  from  the  atmosphere,  density  and  composition  of  the  at,.— spiere,  aurorae 
and  airglow;  magnetoepheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  mag-  rtosphere;  solar  physics,  studies  of  solar  magnetic 
fields;  space  astronomy,  x-.  ly  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia- 
tiona  in  apace  on  apace  ayatema. 
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